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1. INTRODUCTION 

Two-color QCD has been extensively studied over the years both theoretically ^ g || f§] and 
numerically [||, ^, ^ |8j, ^, 10, [11, 12]. Although a lot of progress has been made in uncovering 
important qualitative features of this theory, many interesting quantitative questions remain: 

1. What is the order of the finite temperature chiral transition at zero and non-zero chemical 
potentials? 

2. Can the low energy physics at small T and pL be captured by chiral perturbation theory? An 
answer to this question was attempted in [pQ. 

3. What is the order of the phase transition that occurs when the lattice gets saturated with 
baryons at T = 0? 

4. What is the phase structure in two spatial dimensions since spontaneous symmetry breaking 
is forbidden at finite temperatures. 

The reason for the lack of quantitative progress can be traced to the fact that all previous studies 
have been limited to small lattice sizes and relatively large quark masses; a problem which haunts 
all numerical studies of strongly correlated fermionic systems. 

In this work we try to make progress by considering the strong coupling limit. Although this 
limit has the worst lattice artifacts, it contains some of the essential physics, namely confinement 
and chiral symmetry breaking. On the other hand recent advances in Monte Carlo algorithms allow 



us to study the chiral limit on large lattices with relative ease in the strong coupling limit [J13p . In 
this work we extend these algorithms and apply it to study strong coupling two-color QCD with 
staggered fermions. This theory is especially interesting due to an enhanced U (2) symmetry at 
zero quark mass and baryon chemical potential. It was originally considered in [[14], 15] and was 



recently reviewed in [16]. However, many of the questions raised above remain unanswered even 
in this simplified limit. 



2. THE MODEL 



Our model lives on a d + 1 dimensional hyper-cubic lattice with sites x = (x t ;x\,X2, --^d)- The 
size of the lattice is taken to be L t x L d and is periodic in all directions. The action of our model is 

S = -Y, r a T la( x ) 

x,a 

The Grassmann fields %{x) and %(x) represent row and column vectors with 2 color components 
associated to the lattice site x. The color component of the quark fields will be denoted as % a , a = 
1,2. The gauge fields U a (x) are elements of SU (2) group and live on the links between x and x+ a 
where a = 1 , 2 , . . , d for spatial links and a = t for temporal link. The factor r a = 1 for a = 1 , 2, . . , d 
and ft — ji with a t being the asymmetry factor between spatial and temporal lattice spacing. This 
asymmetry allows us to study finite temperature behavior [|17|| . 



X{x)U a (x)x(x + a)-e-^ 5 <«x(x + a)ui{x)x(x. 



(2.1) 
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A discussion of the relevant symmetries of the action ( |2.1[ ) can be found in [j8|, 16]. As ex- 
plained in these references when /i = our model has a U (2) global symmetry: 

X ^VX , X e ^X e V\ V = exp(ice-d + i0) € 17(2). (2.2) 

where X e and X are given by 

X e = (x e -XT^2), X =( *° ) (2.3) 

and the subscripts e and o refer to even and odd sites. Note in our notation a are Pauli matrices 
that mix % an d %"" present in X a and X e while % are Pauli matrices that act on the color space. The 
U (2) symmetry is reduced to t/s(l) X U x (l) in the presence of a chemical potential: 

t/ B (l): X„ ->exp(iCT 3 0)X o , Z e ->-X e exp(-K730) 
f7 z (l): X ->exp(i0)X o , X e X e exp(-/0). 

Here is the baryon number symmetry %(x) — > e'^^(jc), ^(jc) — ► x( x ) e ~"^ and ^(1) i s 

the chiral symmetry of staggered fermions %(x) — ► Q l ^ e ( x >%(x), %(x) — > x{x) e ' where e(jc) = 
( — i y x i + x i +- v 2 + • • ■ + x <i 

3. DIMER-BARYONLOOP REPRESENTATION 

One of the computational advantages of the strong coupling limit is that in this limit it is 
possible to rewrite the partition function, 

Z = J[DU][dxdx]exp{-S), (3.1) 

as a sum over configurations containing gauge invariant objects [[ll|, 19], 2C]. In our case these 
objects turn out to be dimers and baryonloops. A lattice configuration K of dimers and baryonloops 
is constructed as follows: 

(a) Every link of the lattice connecting the site x with the neighboring site x + a contains either 
a dimer k a (x) = 0, 1 , 2 or a directed baryonic bond b a (x) = — 1 , 0, 1. b a (x) = 1 indicates the 
direction is from x to x + a and — 1 implies it is from x + a to x. k a (x) = and b a (x) = 
means that the link does not contain any dimer or baryonic bond. In our notation we also 
allow a to be negative. Thus, if a was positive, k- a (x) and b- a (x) will represent dimers 
and baryonic bonds connecting x with x — a. 

(b) If a site is connected to a baryonic bond then it must have exactly one incoming baryonic 
bond and one outgoing baryonic bond. Thus baryonic bonds always form self-avoiding bary- 
onloops. 

(c) Every lattice site x that does not contain a baryonic bond must satisfy the constraint 

a 

where the sum includes negative values of a. 
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Figure 1: An example of a dimer-baryonloop configuration. 

An example of a dimer-baryonloop configuration is shown in Figure [[} Note that sites connected 
by k a (x) = 1 also form loops. Given the set {K} of such dimer-baryonloop configurations the 
partition function of the theory described by eq.( |2.l| ) can be rewritten as [|l5|], 



{k} { * 



(3.2) 



where T 



Note that the partition function has been written as a statistical mechanics of dimers 



and baryonloops with positive definite Boltzmann weights. It is possible to extend the Monte Carlo 



algorithm developed in [ |13| ] and apply it to this problem. The details of the algorithm will be 
published elsewhere. 



4. OBSERVABLES 

A variety of observables can be measured with our new algorithm. We will focus on the 
following: 

(a) The chiral two point function, given by 



Gc(z,z') = ( %{z)x(z) J(z')x(z') 



(4.1) 



and the chiral susceptibility, 



(4.2) 



where £2 is the lattice space-time volume. 
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(b) The diquark two point function, given by 

GbM = (xi(z)X2(z) X 2 V)Xi(z')\ (4-3) 
and the diquark susceptibility, 



Xb^^G b (z,z') (4.4) 



(c) Baryon density, defined as 



1 dlnZ 



Y c = 7K £ UE»*M) 2 ) (4-6) 



(d) The helicity modulus associated with the U(l) chiral symmetry, which we define as 

— E 

a=\,2,..,d \ s x 
where 

A„(jc) = e(jc)[|fe (jc)|+*a(jc)] (4.7) 
and Q. s is the spatial lattice volume. 

(e) The helicity modulus associated with the U(l) baryon number symmetry, which we define 
as 




_ (4-8) 
alis a=\X..4 \ x x / / 

where 

B a (x) = [b a (x)} (4.9) 

Both Yc and Yg are diagonal observables and can be calculated configuration by configuration in the 
dimer-baryonloop language and averaged. On the other hand Gc(z,z') and Gb(z,z') are examples 
of off-diagonal observables and can be measured by exploiting the special properties of the directed 
loop update 



5. EXPECTED PHASE DIAGRAM 

At jtt = 0, one expects a finite temperature phase transition separating the low temperature 
(U (2) — ► U (1)) broken phase and the high temperature symmetric phase. Similarly, for small /I 
there is a phase transition separating the low temperature phase where C/b(1) <8>U x (l) is broken 
completely and the high temperature phase which is symmetric. The order of both these transi- 
tions remain unclear. At zero temperature as \i increases one expects the lattice to get saturated 
with baryons which leads to a phase transition from a super-fluid to a normal phase. Renormaliza- 
tion group arguments suggests that if this phase transition is second order it will be a mean field 
transition^]. 

Little is known about the phase diagram in two spatial dimensions. One possible phase struc- 



ture was discussed in [23] in the context of the continuum theory using an effective field theory 
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approach. However, due to infrared divergences that occur at finite temperatures a complete pic- 
ture could not be inferred. Although continuous symmetries cannot break in two spatial dimensions 
at finite temperatures, a Berezinky-Kosterlitz-Thouless (BKT) type phase transition can exist 

6. RESULTS 

Although spontaneous symmetry breaking cannot occur in a finite volume, one can still con- 
clude that the symmetry is broken in the infinite volume limit by studying the behavior of vari- 
ous observables as a function of the volume. In our case the U (2) symmetry at /i = implies 
Gc(z,z') = 2Gb(z,z') and so %c = 2%b- In the presence of a chemical potential since the U(2) 
symmetry is broken, this equality no longer holds. The formation of a diquark condensate can be 
inferred from the growth of %b with the volume. Further, the helicity modulus Y c and Y B , both 
must reach a non-zero constant if the U x (l) and £/g(l) are broken. All these expectations can be 
understood quantitatively using chiral perturbation theory in the £-regime based on an effective 
action, which at jj. = turns out to be 



S e & = I d d x 



where u(x) is a unit two-vector field and S(x) is a unit three vector field. This chiral Lagrangian is 



equivalent to other chiral Lagrangians found in the literature [11]. However, we note that the fact 



that F and B may not be the same was not considered in earlier work . Using the chiral Lagrangian 



is straight forward to extend the results of [ ]25[ ] to obtain a finite size scaling formula for various 
quantities. 

6.1 d = 3 with n = 

The finite temperature phase transition that restores U(2) — > U(l) symmetry breaking can be 
studied in our model by tuning T at fixed L t . We have performed extensive calculations at a fixed 
L f = 4 for different spatial lattice sizes L varying from 16 to 256 and for many different values of 
T. We look for two signatures of the broken phase: 

(a) Both Yc and F# must go to non-zero constants at large L. These constants are equal to the low 
energy constants, F 2 and B 2 in eq.(5.1), of a three dimensional low energy effective theory. 



We use the relations 

Y c = F 2 + b/L + c/L 2 ; Y B = B 2 + b' /L + c /L 2 . (6.2) 
to extrapolate our data to extract F 2 and B 2 . 
(b) The finite size scaling of the chiral susceptibility Xc can b e shown to be 

Zc = y{L 3 + /3 1 (A + J_) L 2 } + aL (6.3) 
where j8i = 0.226 is the shape coefficient for cubic boxes. 
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Figure || shows our results at T = 2.918, a point in the broken phase. As can be seen from the 
graph, the above expectations are satisfied extremely well. In particular we find F 2 ^B 2 . 

Figure [| shows the dependence of %c as a function of L for different temperatures. While 
%c increases as L? at T = 2.9275, it saturates at T = 2.9285 and 2.9292 for large L. Note that 
extremely large lattices are necessary to see this saturation. Thus, there is a phase transition such 
that 2.9275 <T C < 2.9285. The peculiar non-monotonic behavior of %c before it saturates appears 
to rule out a second order behavior but is consistent with a first order transition. We can fit %c to 
the form 

A + BL? exp ( — AFL? ) 
XC = 1 + Cexp(-AFL3) (6 " 4) 

which can be motivated as arising due to the presence of two phases whose free energy densities 
differ by AF. We find this form captures the structure in the data well for L > 48 for both T = 
2.9285 and 2.9292. 

In the high temperature phase one can compute screening lengths E, by looking at the ex- 
ponential decay of Gc(z,z') for large spatial separations between z and z' ■ Similarly in the low 
temperature phase F 2 and B 2 have dimensions of inverse length and hence can be used to provide 
natural length scales in the problem. At the phase transition we find that none of these length scales 
diverge but all are of the order of 40 to 50 lattice units indicating that the transition is a rather weak 
first order transition. 

A renormalization group analysis of the fluctuations of the order parameter field, which in our 
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Figure 3: Plot of Xc versus L for different values of T across the phase transition. The solid lines at 
T = 2.9285 and T — 2.9292 are fits to the first order form discussed in the text. 



case is a complex three vector field <&j(x),z = 1,2,3, has been performed using the £-expansion 



26|, |27|]. Interestingly, the universal field theory one studies also describes the possible normal- 
to-planar super-fluid transition in 3 He flffi]. One finds no stable fixed point which implies that any 
observed phase transition must be a fluctuation driven first order transition. Our observations favor 
this conclusion. On the other hand recently it has been proposed that the £ -expansion results may 
be misleading [|28|]. It is well known that first order transitions can arise due to model dependent 
features. In order to minimize such dependences it may be useful to repeat the above calculation 
for other values of L t . 



6.2 d = 3 at zero temperature 

In order to study zero temperature results we set L t = L, T = 1 and compute observables for 
various values of pL and L. In this case since C/g(l) <8> U x (l) is expected to be broken completely 
the following should be observed: 

(a) The diquark susceptibility %b should grow with the volume, 




(6.5) 



where A = (X1X2) = {%2%\) 7^ is the diquark condensate. 

(b) The chiral susceptibility %c should saturate with L showing that the Gc(z,z') decays expo- 
nentially for large separations between z and z! . Note that the two correlators Gc and Gb are 
no longer related by a symmetry when jj. / 0. 
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Figure 4: The left plot shows Xb an d Xc as a function of L at ji = 0.01. As expected Xb grows as L 4 while 
Xc saturates. For the same value of jx the right plot shows both Yq and Yg grow linearly with L. 



(c) Both the helicity modulus Yq and Fg, as defined in eqs.( |4.6| ),(|4.8|) should grow linearly with 
L. 

These expectations are borne out in our calculations as can be seen in Figure |J 

As the chemical potential increases the lattice is filled with baryons. This can be seen in the 
inset of figure || where the baryon density is plotted as a function of /x. Baryon super-fluidity 
cannot exist beyond saturation and this leads to a phase transition to the non-super-fluid phase. 
Renormalization group arguments show that this phase transition must belong to the mean field 



universality class for d > 2 [|22|]. It was recently shown that the critical chemical potential jj. c = 
0.5 cosh -1 (-/TO) =0.909223.. in the mean field approximation []16[]. The diquark condensate was 
also shown to be 1 

^(\/l0-cosh(2/i)) (6.6) 
In our calculations we extracted A by fitting %b to the relation 

A 2 

Xb = —[L 4 +AL 2 + B}. (6.7) 

Figure |5] shows our results along with the mean field result and the result with one-loop corrections. 
We find that pL c is in excellent agreement with mean field theory while A requires the inclusion of 
one-loop corrections. 



'There is a factor of two mismatch the formula quoted here and what can be found in Jig], The origin of this 
mismatch is the normalization of our kinetic term in eq.([2. l|) as compared to the one in 
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Figure 5: The plot of A and <^ _1 as a function of fi. The dotted line is the mean field result and the solid 
line includes the one loop corrections. The dashed line is a fit to the form =AW([l — fl c ). The inset 
shows the baryon density ng as a function of /!. 

For jti > jU c it costs energy to remove a single baryon and we expect this energy to grow 
as (/» — pic)- Since this phase describes non-relativistic particles the spatial correlation length 
obtained from Gb(z,£), must scale as (jU — jU c ). Figure |5[ also shows that this expectation is 
borne out. 

6.3 d = 2 Phase Diagram 

We have also studied both the zero temperature and finite temperature phase transitions in two 
spatial dimensions. While the zero temperature transition is a mean field transition from a super- 
fluid phase to the saturated phase like in d = 3, the finite temperature phase transition at /» = is 
different and interesting. Mermin- Wagner theorem forbids the breaking of a continuous symmetry 
in d = 2 at finite temperatures. However, the U(l) part(s) of the U (2) symmetry at ju = can 
undergo a BKT type phase transition. This can result in long range correlations in %c an d %b at low 
temperatures. If this is true one would expect 

v \ Const. T < T c 

T>T C (68) 

while lim^co Yb = 0. Further, the screening lengths obtained from Gc(z„z') will behave as 

i=^v(y=), T>T, (6.9) 

Figure ^| shows that indeed our data is consistent with these expectations. We see that Yq is a 
constant and Y B decreases (although very slowly) with L at T = 0.1. For T > 0.7 we can begin to 
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Figure 6: 



see both Yq and Yb decrease with L. For 0.8 < T < 1.03 E, extracted from Gc{z,z. f ) fits well to the 
BKT form with T c ~ 0.6. In a typical BKT phenomena, the super-fluid density is expected to show 
a universal jump at the transition. A naive estimate of this jump suggests that Yq{T = Tc) = 2/n ~ 
0.6366..., in our normalization. Our data shows a different jump suggesting that the transition is 
quantitatively different although qualitatively similar to the familiar BKT transition. 
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